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Résumé. Nous présentons dans cet article une nouvelle architecture optoélectronique pour la génération de chaos 
en intensité. Le principe s'appuie sur une dynamique électro-optique non linéaire à retard, dont la non linéarité est 
construite grâce à un interféromètre à 4 onde réalisé en optique intégrée, et disposant de 2 électrodes de modulation 
indépendantes. Le montage permet de disposer d'une part, d'une dynamique ultra-rapide jusqu'à des fréquences 
de plusieurs GHz, et d'autre part, de générer un chaos de grande dimension destiné au cryptage physique de 
données optiques. Au travers d'une étude numérique et expérimentale, nous avons cherché à analyser certains 
des nombreux comportements dynamiques que peut présenter cet oscillateur, en fonction de divers paramètres 
physiques du montage : régimes de point fixes stables, périodiques, et chaotiques. La mise en œuvre du montage 
expérimental a permis de valider le modèle théorique adopté pour les simulations. 

Abstract. We présent a new optoelectronic architecture intended for chaotic optical intensity génération. The 
principle relies on an electro-optic non-linear delay dynamics, which non linearity is performed by a 4-waves 
integrated optics interferometer involving 2 independent electro-optic modulation inputs. The setup allows both 
to have an ultra- fast dynamics up to several GHz frequencies, and potentially high dimensional chaos intended 
for encryption of optical data at the physical layer. We have built a mathematical model of the System and 
analyzed a number of its possible solutions : stable steady states, periodic and chaotic régimes. The expérimental 
observations allowed to validate the dynamical model, through good qualitative agreements with the numerical 
simulations. 

1 Introduction 

Les communications sécurisées par chaos ont permis le développement de nombreux systèmes dyna- 
miques optiques produisant des comportements chaotiques complexes [1]. La plupart de ces réalisations 
utilisent une classe particulière de comportements chaotiques, celle des dynamiques non linéaires à retard 
[2]. Ces dynamiques temporelles possèdent la particularité étrange d'évoluer dans un espace des phase de 
dimension infinie, dans lequel des comportements chaotiques de grande complexité peuvent être observés 
[3]. L'oscillateur électro-optique que nous allons présenter fait partie de cette catégorie, et le chaos est 
observée sur la variable intensité optique. 

Tout d'abord, le dispositif expérimental sera décrit, ce qui permettra d'en établir un modèle théorique 
sous forme d'un système de deux équations différentielles couplées du second ordre, non linéaires, et à 
retards multiples. Ce modèle sera par la suite utilisé pour effectuer des simulations numériques. Ensuite, 
nous présenterons et comparerons les résultats obtenus par des mesures expérimentales et ceux obtenus 
numériquement. Enfin, nous conclurons en décrivant le système cryptographique complet, qui doit à terme 
utiliser ce générateur de chaos comme dispositif émetteur d'un côté, et comme récepteur / décodeur de 
l'autre côté. 



2 Générateur de chaos à modulateur QPSK 

2.1 Description et principe de fonctionnement 

Le dispositif expérimental de l'oscillateur chaotique est illustré sur la figure 1. Le générateur est formé 
par 2 boucles de contre-réaction reliées à un modulateur QPSK (Quadrature Phase Shift Keying). Ce 
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modulateur électro-optique est particulier par son architecture. Il appartient à la famille des modulateurs 
Mach-Zehnder (MZ) à 4 bras, et il est intégré sur Niobate de Lithium (LiNbOs). Cet élément clé est un 
composant commercial récent, originellement destiné à des nouveaux formats de modulation numérique 
pour les télécommunications optiques [4]. Il permet pratiquement dans notre cas la réalisation d'une 
fonction non linéaire bidimensionnelle. 
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Fig.l. Schéma du dispositif expérimental d'un générateur de chaos à modulateur QPSK. 



Le reste du générateur de chaos est plus classique. Il est constitué (de gauche à droite) par les composants 
suivants : 

- une diode laser monomode reliée à une fibre optique à maintien de polarisation ; sa longueur d'onde 
d'émission est celle habituellement utilisée dans les télécommunications optiques (Ao = 1,55 fim). 

- le modulateur électro-optique QPSK déjà éviqué. Il a la particularité d'avoir 2 entrées RF indépendantes 
et 3 tensions de polarisation indépendantes. Ce composant représente le cœur du système, sa fonc- 
tion de transfert de modulation réalisant la fonction non linéaire du générateur de chaos. 

- le système présente aussi l'originalité d'avoir 2 boucles de contre-réactions électro-optiques non 
linéaires, impliquant chacune un élément retardant, une photodiode de conversion optique / électrique, 
un filtre électronique large bande responsable des termes différentiels de la loi, et un élément am- 
plificateur. 

1. La première boucle est réalisée par la mise en série des éléments suivants : 

- un coupleur optique (2 entrées / 2 sorties) permettant d'insérer le message utile par l'une de 
ces entrées, et d'envoyer le mélange (chaos + message) d'une part vers la première contre- 
réaction, et d'autre part vers la canal de transmission. 

- une fibre optique réalisant un retard temporel pur, et constituant un paramètre physique 
« clé » de la première boucle de contre-réaction (La valeur très précise de ce retard est 
nécessaire pour l'opération de synchronisation / décodage). 

- une photodiode permettant de convertir la puissance optique en un signal électrique. 

- un filtre passe bande limitant la bande passante du signal chaotique électrique. 

- un amplificateur électronique réglable pour ajuster le gain de boucle de la contre-réaction. 

2. La deuxième boucle de rétroaction est réalisée par des éléments semblables à ceux utilisés 

dans la première boucle, à l'exception du coupleur optique à (2 entrées / 2 sorties) qui ne fait 
plus partie de cette seconde boucle de contre-réaction. 

Le processus dynamique réalisé par cet oscillateur en double boucle fermée peut être décrit de la 
manière suivante : une diode laser fibrée alimente optiquement le modulateur QPSK. Le faisceau lumi- 
neux traversant le modulateur QPSK subit à la sortie une variation d'intensité non linéaire par rapport 
à chacune des deux tensions appliquées aux électrodes de commandes RFl et RF2 du modulateur. Le 
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coupleur (1 entrée/2 sorties) divise ensuite cette intensité optique en deux quantités égales. Les signaux 
issus des sorties de ce coupleur sont ensuite retardés différemment dans chaque boucle par une certaine 
longueur de fibre optique (quelques centimètre à quelques mètres), puis ils sont détectés par des photo- 
diodes. Les signaux électriques obtenus sont filtrés, puis amplifiés pour y être réinjectés sur les électrodes 
de modulation RF du QPSK. La présence du coupleur optique (2 entrées / 2 sorties) dans la première 
boucle ne sert qu'à mélanger le message informatif avec la modulation chaotique de l'intensité. 

2.2 Modélisation de la fonction non linéaire 

La figure 2. a schématise la configuration du modulateur QPSK utilisé pour la réalisation de la fonc- 
tion non linéaire bidimensionnelle, notée par : fjsfL \^a-> Vb\ • Cette dernière est représentée par la fonction 
de transfert du modulateur. L'expression de l'intensité optique, qui représente la variable de sortie du 
modulateur QPSK, est fonction des tensions appliquées sur ces électrodes RF (v a sur RF1 et sur RF2) 
ainsi que de trois tensions continues (DCl, DC2 et DC3) qui permettent d'ajuster des points de repos de 
la condition d'interférence à onde multiple. 




(a) Schéma élémentaire (b) Modélisation (c) Allure de la fonction de transfert 



Fig.2. La fonction non linéaire réalisée par le modulateur QPSK. 



L'expression analytique de la fonction de transfert du modulateur QPSK peut s'établir de la manière 
suivante : 

- Sur l'électrode DC m (m = 1,2,3) est appliquée une tension continue Vdc^- Un déphasage (j) m est 
introduit sur l'onde optique qui se propage le long de la branche soumise à la tension Vbc m (comme 
illustré sur la figure 2.b). L'expression de ce déphasage est donnée par : 

Vnc 

m = 7T.-^- (1) 

où V^DCm est la tension demi-onde, qui permet de réaliser un déphasage de tt. 

- Sur les électrodes RF12 sont appliquées respectivement 2 tensions variables v a ,b(t)- Ces dernières 
sont des tensions de modulation, et elles introduisent des déphasages variables ipi^ (t). Les expres- 
sions de ces déphasages sont données par : 

<Pl,2{t)=^-^- + 'h,2 (2) 
VirRFi t 2 

où V 7V rf 1 2 sont des tensions demi-ondes, qui permettent de réaliser en régimes dynamiques un 
déphasage de tt respectivement dans les interféromètres de Mach-Zehnder MZi et MZ2. 
L'expression de l'intensité optique I ou t(t) à la sortie du modulateur QPSK est obtenue par le calcul de 
la moyenne (sur un temps défini par le temps de réponse électronique) du module au carré du champ 
électrique résultant E out {t) : 

Iout(t) = (\E out (t)\ 2 ) = f NL [v a ,Vb](t) (3) 
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En imposant un champ en entrée du modulateur sous la forme Ei n (t) — \/ r Pô. ex.-p(icoot) (Po est la 
puissance optique à l'entrée du QPSK, c^o = 2tï.c/\q est la fréquence angulaire de la source laser, et c est 
la vitesse de la lumière), le champ électrique E out (t) en sortie du modulateur QPSK est donné par : 



E ou t(t) 



fPÔ 



1 + exp (iipi(t)) + 1 + exp (i(p 2 {t)) • exp (i0 3 ) 



exp (iwot) 



(4) 



Il vient finalement pour la fonction de modulation non linéaire bidimensionnelle 



fNL[v a ,Vb](t) = y<j C0S(^3) ' 



COs('03)+2-COs('0i+'02)-COs('02+'03-'0l) + COS 2 (^2 + ^3 " ^1 ) | (5) 



avec : 
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^2 
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^3 



Une condition suffisante sur la non linéarité pour obtenir une dynamique chaotique est de présenter un 
extrémum dans l'intervalle de variation des variables d'entrées. Cette condition est suffisamment vérifiée 
par fjsfL \pai Vb\ i donnée par la relation (5), comme le montre un exemple de son allure sur la figure 2.c 
(les paramètres utilisés pour tracer cette figure sont donnés au tableau 1). 



3 Modélisation du système 

Le processus différentiel se modélise à partir de la fonction de transfert du filtre passe-bande de la 
branche électronique de l'oscillateur. Comme l'oscillateur est formé par 2 boucles de contre-réactions 
contenant 2 filtres passe-bande différents, nous adoptons pour désigner les paramètres de chaque boucle 
les notations suivantes : tous les paramètres de la boucle (A) seront indexés par la lettre « a » , et tous les 
paramètres de la boucle (B) seront indexés par la lettre « b ». Ainsi, la dynamique globale est modélisée 
par 2 équations différentielles du second ordre à retard (système d'EDR), dont leurs expressions sont 
données par : 



ÎNL[x a ,Xb}{t - Ti) 



(6) 



où (i = a, b) selon la boucle concernée, et xi(t) = Vi(t)/(2 • V^rfï) représentent les variables normalisées 
où Vi(t). r\i et T2i sont les constantes de temps caractéristiques du profil des filtres électroniques, elles 
sont reliées directement aux fréquences de coupures f c u et f C 2i haute et basse du filtre, = (tt ■ Po ■ 
7o • Gi • Si • ao • cti)/(2 • V^rfï) est le gain global normalisé de la boucle de rétroaction, avec : 70 le 
coefficient des pertes optiques du modulateur QPSK ; Gi le gain de l'amplificateur ; Si la sensibilité du 
photodétecteur; ao le coefficient de couplage du couleur optique (1x2). est le coefficient de couplage 
du coupleur optique (2 x 2) (qui n'apparaît pas pour la boucle (B)). 



4 Résultats numériques et expérimentaux 

Les simulations ont été effectuées en intégrant le modèle donné par (6) par la méthode du prédicteur- 
correcteur. Un outil important dans la compréhension de la route vers le chaos est le diagramme de 
bifurcation, qui permet de voir rapidement l'ensemble des régimes dynamiques obtenus pour différentes 
valeurs du paramètre de bifurcation. Ce dernier correspond dans notre cas au gain normalisé (3 a de la 
contre-réaction de la boucle (A). 

Sur les figures 3. a et 3.b, nous avons représenté les diagrammes de bifurcation obtenus expérimentalement 
et par simulation numérique. Afin de faciliter la comparaison expérience / simulation, le mode de 
représentation correspond à l'évolution, en fonction d'un gain de boucle, de la densité de probabilité 
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de la variable normalisée x a (résultat du filtrage passe-bande de fjsfL [^a? %b] •> dans la boucle (A)). 
Des premiers diagrammes ont été effectués dans le cas d'une seule boucle de rétroaction, cas qui peut 
servir de référence puisqu'il est déjà étudié dans la littérature [5]. Le système en une seule boucle signifie 
que physiquement, une boucle est laissée ouverte, et théoriquement, il se traduit par un gain nul (/?& = 0). 
On remarque sur ces diagrammes l'évolution typique de la dynamique du système jusqu'au chaos selon 
une route du type cascade par dédoublements. Les dynamiques chaotiques obtenues se traduisent par une 
entropie élevée, et une densité de probabilité diffuse, et de profil quasi-Gaussien. 



Boucle (A) 


Boucle (B) 


symbole valeur unité 


symbole valeur unité 


T a 61 ns 
fcla 13 GHz 

fc2a 50 kHz 
Tla = 2.Jf cla 12 > 2 P S 
T2a = 2.7r./ c2o 3 ' 18 ^ S 


T b 60 ns 
fdb 13 GHz 
Ub 30 kHz 

™ = 2. W ./ Clb 12 > 2 P S 
T 2b = 2 .K.f c2h 5 > 30 M S 


Paramètres de la non linéarité 


V*RF1 5,84 V 
V^Dci 7,40 V 
V*dc2 7,14 V 
V*DC3 14,24 V 


V*rf2 6,08 V 

01 2,9 rad 

02 1,3 rad 

03 -0,1 rad 



Tab.l. Paramètres expérimentaux du générateur de chaos à modulateur QPSK. 




(a) expérimental (b) simulation (c) entropie 



Fig.3. Diagrammes de bifurcation et entropique du système en une seule boucle de rétroaction (/3b = 0). 




(a) expérimental (b) simulation (c) entropie 

Fig.4. Diagrammes de bifurcation et entropique du système en doubles boucles de rétroaction = 0, 1). 
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Les diagrammes de bifurcation du système en double boucle sont donnés sur les figures 4. a et 4.b. 
On remarque sur le diagramme de bifurcation expérimental (figure 4. a) les différents régimes dynamiques 
(régimes points fixes stables, régimes périodiques et chaotiques). Le diagramme obtenu par simulation 
reproduit le résultat expérimental avec une très bonne ressemblance (figure 4.b). 



5 Conclusion 

Dans cet article, nous avons présenté une nouvelle architecture de générateur de chaos en intensité, 
basé sur une dynamique non linéaire à retard. Le dispositif expérimental est réalisé par des composants 
électro-optiques, dont le cœur du système est un modulateur QPSK. Ce modulateur permet de relier 2 
boucles de rétroactions, et d'introduire une non linéarité en 2 dimension. Cette architecture du générateur 
nous permet de disposer d'un nombre important de paramètres physiques (clés de codages) pour réaliser 
des communications sécurisées au niveau physique à haut débits. 

Nous avons modélisé le système par un système d'équations différentielles du second ordre à retard. 
La comparaison des résultats de l'intégration de ce système par la méthode prédicteur-correcteur aux 
résultats expérimentaux, nous a permis de valider le modèle théorique proposé. Dans les simulations, 
nous avons pris soin de choisir un pas d'échantillonnage constant (h = 1 ps) assez petit devant la plus 
petite constante de temps du système {r\ a = T15 = 12,2 ps). Nous avons choisi aussi pour la réalisation 
du générateur, des retards temporels très grand devant la plus petite constante de temps afin de garantir 
la génération d'un hyperchaos (T^/rib ~ 5000 ^> 1), et d'optimiser par conséquent la confidentialité des 
transmissions. 

Le travail va désormais se poursuivre sur la réalisation du système global de cryptographie par chaos, 
dont les premiers résultats numériques sont déjà obtenus en modulant chaotiquement une information 
binaire à plus de 3 Gbit/s. La restitution du message au niveau du récepteur est basée sur le principe de 
la synchronisation de chaos, initialement introduit par Pecora et Carroll [6]. 
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